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Abstract 

It has been conjectured0'i that an impurity with charge Z >2 can be locahzed 
due to its interaction with electrons in a metal. The simplest case is an 
impurity free to move between only two sites, which interacts locally with 
s-wave electrons. For Z >2 the hopping of the impurity is formally irrelevant 
and this has been argued to lead to localization. In this paper it is shown that 
other processes, in particular joint hopping of the impurity and one or more 
electrons between the sites, play an important role and have not been treated 
properly in the literature. Being relevant in a renormalization group sense, 
even when Z > 2, these terms lead to delocalization of the impurity. Using 
bosonization, it is shown how these processes are generated from marginal 
operators that are usually neglected and the dangers of ignoring marginal or 
irrelevant operators are discussed in detail. Questions about implications for 
the more general situation of many sites to which the impurity can hop, are 
also considered. 
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I. INTRODUCTION 



A very interesting but subtle problem for the past quarter century has been the low 
temperature (T) behavior of a heavy particle interacting with the electronic excitations in a 
metal. The subtleties arise from the competition between the tunneling of the heavy particle, 
which tends to delocalize it, and the inability of the electronic degrees of freedom to adjust to 
the potential of the moving particle, which tends to localize it. This difficulty is manifested 
in the orthogonality between the electronic ground states with the impurity located at two 
different points in space. This phenomenon, Anderson's orthogonality catastropheBi, is due 
to the fact that the two ground states differ by a very large number of very small energy 
particle-hole excitations. This infinity of excitations is related to the fact that the impurity 
at two different points in space creates Friedel oscillations in the electronic charge density of 
the metal which, due to a difference in the phase of the oscillations, differ from each other 
at arbitrary distances, implying particle-hole deformations at arbitrarily low energies. 

In the limit of weak interactions between the heavy particle and the electrons the primary 
effect of the coupling to the electrons is to induce a frictional force on the otherwise free 
tunneling motion of the particle — although the dynamic properties of the impurity at low 
temperatures are only partially understood!. In the opposite limit of strong interactions, it 
has been argued that the particle will be strictly localized. More specifically, it has been 
claimed that a charged particle with charge Ze > 2e tunneling between symmetric positions 
that are well separated spatially will localize around one of the sites due to the interaction 
of the particle with s-wave electrons. 

In this paper we will argue that a heavy particle interacting with the electrons via a small 
number of channels (less than or equal to four) cannot be localized by the interaction because 
of subtle effects that have not been treated fully in previous workllli'iS. We will primarily 
focus on the simplest case in which the particle can tunnel between only two sites that are 
related by symmetry. This two site problem is related to the Kondo problem although there 
are important differences which have often been ignored in the literature. We primarily 
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assume local, screened, s-wave, spin independent interactions between the particle and the 
electronic degrees of freedom. Due to the spin independence of the interaction we can then 
neglect any possible spin of the heavy particle and treat the opposite spin electron species 
independently. 

At the end of the paper we will discuss the generalization of the problem due to the 
presence of three or more sites to which the particle can hop, and also the potential relevance 
of more angular momentum channels. We will argue that our results suggest the correct 
behavior for the extended system of a particle in a periodic metal, and also have implications 
for the sharpness of X-ray edge singularities in systems with mobile deep holes and other 
related problems. The purpose of this paper is partially pedagogical, thus we work through 
some parts in substantial detail, in particular pointing out the dangers that lurk within 
many of the standard tricks, especially bosonization. 

A. Outline 

In the remainder of this Introduction we motivate the form of the Hamiltonian with which 
we will primarily work and explain qualitatively the effects of the orthogonality catastrophe 
on the motion of the particle as well as the effects that make it difficult to localize. The rest 
of the paper is organized as follows: In the next section (II) we motivate and introduce the 
standard and very useful method to perform the calculations, i.e. bosonization. The model 
is introduced in the usual fermion representation of electrons which is then mapped into 
bosons. In section III a path integral representation of the partition function is formulated 
and brought into a Coulomb gas representation. Renormalization group flow equations 
are derived and analyzed. In this way the results that were discussed qualitatively in the 
Introduction are put on a firm footing. Subsequently in Section IV, the results, as well 
as possible generalizations and complications, are discussed. Finally in Appendix A the 
Coulomb gas representation of the partition function is rederived from the original fermion 
representation and in Appendix B the two site problem is analyzed in the absence of any 
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symmetries other than the equivalence of the two sites. 



B. Physical Picture 

The Hamiltonian of the impurity-electron system has three important terms: the non- 
interacting electron part {Ho), the hopping of the heavy particle between the two sites and 
the interaction term U. Thus: 

H = Ho + A„(ci+4 + d+di) + U (1) 

where df,d2 are the creation operators of the impurity at sites 1, 2 and Ag is the bare 
hopping matrix element of the impurity between the sites. We will generally neglect any 
asymmetry between the two sites. The interaction U will involve terms of the form didic^Cn' 
and d2d2C^Cn' where are the creation operators of the electronic degrees of freedom. 
However, due to the assumed local nature of the potential we can rediagonalize the degrees 
of freedom of the electrons and be left with only two electronic degrees of freedom for each 
energy that are just composed of those wavefunctions which do not vanish at the two sites. 
Thus the potential U can be put in a form in which it involves two electronic states — albeit 
not free electron eigenstates — while all others decouple from the impurity. We can thus 
write the most general form for U in the following symmetric way: 

U = (dtdi + d^d2)[Vi{4ci + c+ca) + ^2(0^02 + 4ci)] 

+ - d+d2) [V^ictci - C+C2) + tVi{ctc2 - C+Ci)] (2) 

with matrix elements Vi to be determined. The interchange symmetry 1 ^ 2 is manifest 
in U. We have picked a basis for the electrons q , i = 1, 2 so that in the limit that R, the 
distance between the two sites, tends to infinity, the Cj's tend to the local s-wave annihilation 
operators at the two sites. As a result we expect V2 ,V4 — * and Vi as R 00. This 

will be seen explicitly later. Because there is only one heavy particle it is convenient to 
express U in form of Eq.(^) since dfdi + d2d2 = 1. Therefore only V3 and V4 couple the 
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impurity to the electrons. One can see that the four terms in Eq.(|]) are the only ones 
possible, due to the 1^2 interchange symmetry of U. From these four terms, by choosing 
the basis q appropriately, one can make one term vanish since there are many ways we can 
choose normalized states that all tend to the local wavefunctions of the two sites as i?, ^ oo. 

The freedom of choice of states is related to a gauge symmetry. If the system is time 
reversal invariant, to which we primarily restrict consideration, then a gauge can be chosen 
to make the Hamiltonian real and hence eliminate V4. Note that, more generally, even in the 
absence of time reversal invariance V4 could be eliminated formally at this point. But other 
operators would appear in the more detailed analysis which cannot be eliminated. Although 
we will not analyze these operators in detail, we will argue why they will not affect our main 
results. 

Thus H can be written as: 

H = Ho + Ao(rf+d2 + d^di) + Viictci + 4C2) 
+V2icic2 + c+ci) + V^id+di - rf+rf2)(c+ci - C+C2) (3) 

using dfdi + (ij(i2 = 1. In the next section we will explicitly derive this form and evaluate 
the Vi's. In the standard manner, one can treat the relevant electronic degrees of freedom 
that comprise Ci and C2 as essentially one dimensional with the magnitude of k playing the 
role of a one dimensional wavevector. 

Let us now try to understand the effects of the potential on the motion of the particle. 
To start, we consider the simple limit with R large so that the V2 term which couples the 
two channels vanishes and Vi = V3. Thus we are left with two independent channels, ci and 
C2, that interact locally with the impurity. Channel 1 electrons interact with the impurity 
when it is on site 1, {dfdi = 1) with interaction strength Vi and do not interact when the 
impurity is on site 2, ((i^(i2 = 1). The opposite holds for channel 2. The important physics 
arises from the effect of the impurity at a given site on the electrons. If, for example, the 
potential is attractive, then the impurity tends to attract electrons towards it in order to 
screen its presence. Thus when the impurity is on site 1, it will tend to shift electrons of 
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channel 1 towards site 1. A rather naive picture of this shifting is the induction of charge in 
a metal close to a positively charged object in order to screen the electric field in the bulk 
of the metal. The induced charge comes from the outer boundaries of the metal and thus 
from arbitrarily far away. 

A better interpretation is in terms of the wavefunctions of the electrons. Due to the 
existence of the attractive potential the wavefunctions far away from the potential center 
look just like the non-interacting ones except for a phase shift. This implies that some extra 
charge density has moved in from far away to screen the impurity. Indeed Priedel's sum rule 
relates the phase shifts 5(, at the Fermi level to the charge, Ze, that is needed to screen a 
charged impurity, 

Z = 2^(2£+l) (4) 

where the sum is over the angular momenta channels £; the sum over the spins yields the 
factor of 2. Thus we can interpret 

n(, = Sij-K (5) 

as the number of electrons per channel that need to be shifted close to the given site in order 
to screen the potential. For the present discussion we will assume that only one angular 
momentum channel £ = plays a role. 

Now since the impurity can move from site to site, in order to understand the effect of the 
interaction on the dynamics of the impurity we need to know the time dependent amplitude 
of a process in which the impurity hops away from a given site for a certain amount of time 
t before it hops back to the previous site. When t is long enough we can view this process in 
the following simple way: until time t — Q the impurity has been at, say, site 1. At i = the 
particle tunnels to site 2 where it remains for a time t before tunnelling back to 1. When 
the particle hops away from site 1, there are n — 8/t: extra electrons of each spin within a 
screening distance from site 1 that will move away as the system evolves to its new ground 
state. Similarly there are n extra holes of each spin near site 2. With the sites far apart, 
the evolution of the s-wave electrons around each site are essentially independent. 
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Following Schotte and SchotteQ we can get a semi-quantitative understanding of the 
amplitude of the hopping process by considering bj'K = n with n an integer. Between time 
zero and t, the n extra s-wave spin- up electrons in channel 1 propagate as in absence of the 
potential that earlier kept them near site 1. To estimate the amplitude of the total process, 
we need to find the matrix element between this evolved state at time t and the ground 
state with the particle back at site 1, i.e. the initial state at time zero. Since the radial 
distance from site 1 of the s-wave electrons can be treated as essentially a one dimensional 
coordinate, we can obtain the t-dependence of this process, roughly, by creating the n 
electrons at distances a, 2a,. . . , raa from site 2 with na of the order of the screening distance. 
(A better approximation would involve an integral over the positions of the extra electrons 
with a weighting factor related to the wave functions in the presence of the impurity at site 
1; but this will only modify our crude estimate by a multiplicative prefactor.) We thus need 
to compute the amplitudelli 

A^{t) = (0|c(a,t) ■ ■ ■ c{{n-l)a,t)c{na,t)c'^{na,0) ■ ■ ■ c+(2a, 0)c+(a, 0)|0) (6) 

in a one dimensional system with no potential, with all the electrons moving at the Fermi 
velocity vp in the same direction. At long times, the antisymmetry of Eq.(^ under ex- 
change of any two space variables fixes the form of A„. When t ^ Tc = a/vp , the sum 
over all possible Wick pairings in Eq.(|]) with the one dimensional long time propagator 
Go oc [i(t — x/vf)]~^ yields 

An ^ det[t + {j ~ t)T4-\ (7) 

The determinant is of the n x n matrix with i and j subscripts. By use of the properties of 
determinants, this can be shown to yield: 

for long t <^ Tc. The same result will obtain for the down spin electrons as well as for 
the s-wave holes around site 2. Thus the amplitude for the full double-hop process will be 



A (xt . In general, with different spin, angular momentum and site channels, 7, with 
electrons moved in channel 7, the amplitude will he A oc t ^ . Later we will see that 
the general result for far away sites is to simply replace by an appropriate phase shift 
= . 

It is convenient to define an exponent ao 

"o = ^E^? (9) 
^ 7 

so that the amplitude of the double-hop process will be, including dependence on the bare 
hopping amplitude Aq , 

~ A^r^"". (10) 

If the sites are not far apart, or the system is not rotationally invariant, there will nevertheless 
still be quantities analogous to n^, with the interpretation as charge moved in a "channel", 
such that Eq.([TO|) obtains, even though the phase shifts no longer have any meaning, see 
Appendix B. 

In order to understand the djTiamics in the presence of the coupling to the electrons, we 
make the standard argument, with the Ansatz that in equilibrium, the heavy particle hops 
back and forth at a rate A.&i The amplitude for this hopping can thus be guessed to be the 
square root of the double-hop amplitude A{t ~ 1/A) since the particle will spend time of 
order 1/A at each site before hopping back. Thus the amplitude will be of order A{t ~ 1/A) 
for each pair of hops, so that, from Eq.(|l^) we have 

A^-A^A^"". (11) 

This has the following solution: 

{0 for > 1 

^ (12) 
cAo'~"° for ao<l 

We thus see that for ao > 1 the real hopping process will not take place and we are thereby 
lead to the conclusion that the impurity will localize on the site on which it started under- 
going only short virtual hops back and forth to the other site. With only s-wave scattering 
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off tfie impurity Friedel's sum rule yields Z = Iboji^ with bo the s-wave channel phase shift 
for the potential ^ = Vi = V3 and the factor of 2 coming from the two spin species. As a 
result, for well separated sites, 



obtaining the inequality if Z > 2, so that a charge two particle will be localized although a 
charge one particle will not be. 

This is the conclusion that has been reached, by this argument and more sophisticated 
versions of it, by a number of author sBi'iil. It seems widely accepted - along with the 
extension of the result to the localization of a particle moving on a lattice of sites. Note, 
furthermore that if there were more angular momentum channels present with phase shifts of 
both signs, it should be possible, by the above argument, to localize even a neutral or charge 
one particle provided the phase shifts are in the regime in which the exponent > 1- The 
main point of this paper is that these conclusions are not justified. Although we will see that 
it still appears to be possible to localize a particle, this cannot be achieved by just s-wave 
scattering for any charge, and in fact requires at least three angular momentum channels to 
have substantial coupling (so that s and the three p channels may be sufficient). 

We will see that the approximation of neglecting the V2 coupling is very dangerous. In 
contrast, relaxing the approximation of V3 = Vi will not change much and the relevant phase 
shifts will be those associated with V3. However the crucial V2 term changes the symmetry of 
the problem: The Hamiltonian in Eq.@ with V2 = and the c\ and ci electrons uncoupled, 
has the continuous extra gauge symmetry 



that is broken by V2. The V2 term mixes the two channels around the two sites (although 
the mixing will be weak for large intersite separations). This term allows processes in which 
one or more electrons near one site transform to electrons near the other. We shall see that 
these yield processes in which the impurity hops from one site to the other simultaneously 




(13) 



(14) 
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with a number of electrons moving from one site to the other. In terms of the interpretation 
of the exponent for the time dependence of a process as a square of the charge transferred 
(such as Eq(^) we see that the exponent for a process in which the impurity and a hole of 
each spin hop together will be 

ai=2K-l)2. (15) 

This process is illustrated in Figure 1. In general the process in which the impurity moves 
from site 1 to 2 at the same time as m holes of each spin transfer from site 1 to site 2, will 
have an orthogonality exponent 

= 2{no - mf. (16) 

Thus using the self-consistent argument for localization outlined above, we conclude that in 
the case of s-wave scattering with no = Z/2, regardless of the charge of the impurity, it will 
never become localized because there will always be a process with m pairs of holes with m 
such that 



f -m 



< I for which 



am = 2 (n„ - mf = 2{-~mf< -. (17) 

This process will yield a non-zero hopping rate and will delocalize the particle. 

Physically, a process in which a number of electrons hop as well as the impurity, schemat- 
ically shown in Fig.l, means that in a sense, less of the screening cloud hops back and forth 
that one would expect from the behavior of the ground state of the static impurity. The 
combined process can be thought of as the tunnelling back and forth not between the static- 
impurity ground states, but between excited states, with the extra action associated with 
this combined process more than compensated for by its larger matrix element (since it has 
a smaller orthogonality exponent). The process with the least action overall will dominate 
the impurity hopping. 

We shall see that this effect can easily be missed, and indeed it seems to have been missed 
in the literatureBSi, even though a number of authorji0 have considered "electron assisted 
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tunnelling" processes in which the heavy particle hops simultaneously with one electron 
(another process that can occur). This process was also introduced in a spinless version of 
the X-ray edge problem in a recent numerical work by Libero and Oliveira0. The main 
theoretical difficulty is that in certain representations (e.g. choices of fields to bosonize) the 
important extra impurity-electron hopping terms are generated, under renormalization, from 
marginal terms (such as V2) which are themselves only generated from irrelevant operators. 
As happens all too often, irrelevant operators cannot just be cavalierly thrown away! 



II. MODEL 

In this Section we introduce a simple model with short-range interactions, show how it 
can be cast in the form of Eq(^ and then begin to analyze it by bosonization of the electron 
fields, pointing out some of the pitfalls. 



A. Definitions 

We start with the Hamiltonian 

H = Ho+Y.Ua + K{dtd2 + dtd{) (18) 

(T=± 

where the free electron Hamiltonian Ho can be written as 

Ho= ^kCkaCka (19) 

with /k = / and with c^o- being the creation operators of electrons at momentum k, 
spin a and energy £k- Finally Ua is a short-range interaction between electrons and the 
impurity 

U, = V [ [ e-'^'^-'^'^-'-ctc^.^dtdi + V I I e-^('^-i''>^^c+,Ckv4^^2 (20) 

Jk Jk' Jk v/k' 

with V the interaction strength and ri the position of the i-th site. If we put the center of 
coordinates between the two sites then we can set ri = — y and r2 = y with R = ri — r2. 
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In order to eliminate the unimportant degrees of freedom that are decoupled from the 
impurity we integrate over the k-solid angleslll and are left with a set of effectively one- 
dimensional degrees of freedom. Neglecting for now the spin index o", we define c^^ via: 



dk f d k j^ik-— -L- ^ dk 



(27r)' 



(21) 



2n J (27r)" " J 27r 

with dQ]f_ being the solid angle element in k-space. But now these one-dimensional Fermi 
operators are not properly orthogonal. This is manifested by nonvanishing anticommutation 
relations {{c^i.,c^k} 7^ 0). An orthonormal set of states can be made from these that are 
even and odd under the interchange of the two sites: 

ctk = ^(Kk - ctk) (22) 



/No 

where the subscripts e,o denote, respectively even and odd while the normalization constants 



e,o 



.r 2A;2 / , sinA;i?\ 

are picked so that Cek and Cok satisfy one-dimensional anticommutation relations: 

{ctk,Cek'} = 27r5ik-k')etc. (24) 
From these states we can obtain linear combinations 

+ _ ^tk ~l~ ^ofc 

4k = (25) 

which transform into each other under interchange of the two sites. It is interesting to note 
that the only orthonormal states that have this symmetry for arbitrary kR. To see 

this one could basically define the most general pair of orthonormal states with interchange 
symmetry: 
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4k = Pe^'cXk + (26) 



sin kR 



Thus for small 



with /c-dependent 60. Eq (p6D is well defined only if |cos6'| > 
kR we must set ^ ~ while for large kR the two sites are decoupled and 6 can take almost 
any value; for 6' = Eq(|26|) becomes Eq(^) with the use of Eq(|22D. Now we can invert 
Eq(l^) and using Eq(|21|), substitute into the potential U in Eq(pOD. Since the only values of 
k that play a significant role are k ^ kp "we can set a, (3 (equivalently Ne, No) to a constant 
evaluated ai kp- Thus using 

4 = 1^0^ with ^ = 1,2 (27) 

with the integral running over k in the neighborhood of kp with an appropriate cutoff of 
order kp, we get an expression for the potential identical to Eq(^. Henceforth we will choose 
^ = for all k which yields V4 = in Eq(|^) thereby explicitly exhibiting the time-reversal 
invariance. 

Furthermore we can obtain the other coefficients Vi in Eq(|^) starting from Eq(^) by 
using the relation between the free electron density of states per spin ai e = ep, pp, with 
the Fermi momentum kp] after rescaling the Fermi velocity to be one, 2tt'^pp = kp. We then 
get: 

Vi = nppV 

V2 = TcppVSi^ (28) 



^3 = r^PFV^l - 

It is also instructive to write U using the even-odd states from Eq(p^). Defining 



we obtain: 

U = Vi (c+Ce + C+Co) + V2 (c+Ce - C+Co) + V3 (c/^rfi - rf^c/a) (c+Co + C+Ce) (30) 
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Comparing Eq(pOD with Eq(^ we see that (with V4 = 0) the non-interacting part of U is 
diagonal in the even-odd representation while the interacting part is diagonal in the Ci^2 
representation. 

Finally, we make the standard change of variables for a two state system, i.e. 

did2 + c^i = cFx- (31) 

In this representation the impurity in site 1 (2) is in state — (+) of the az operator. Thus 
= l~) while (i^ |0) = |+) where |0) is the ground state of Ho- 
Summarizing, the Hamiltonian can be written as 

H = Ho + Ao(T^ + U (32) 

with 

U = Vi (c+ci + c+ca) + V2 (c+C2 + c+ci) + V3 (c+ca - c+ci) (33) 
and with Vi given by Eq(p8D and q given by Eq(^) (where i = 1,2) while 

/dk 
^^kcfkCik (34) 

and 

with Fermi velocity set equal to unity. 

The electrons that interact with the impurity are thus effectively two species of one- 
dimensional fermions moving only to the right, with those to the left of the origin cor- 
responding to incoming electrons while those to the right of the origin corresponding to 
outgoing electrons. 

At this point it is useful to pause and consider the symmetries of the effective Hamiltonian 
in Eq(^). There is a global U{1) gauge symmetry — of the electron phase — and a discrete 
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interchange symmetry 1^2. However note that in the absence of the V2, there would be 
an extra gauge symmetry, that of Eq(|T^. Although in some formulationslE it appears that 
V2 can be made to disappear, this is potentially dangerous as V2 breaks the artificial extra 
gauge symmetry and the formally irrelevant operators which break the symmetry should 
thus be retained. 

It is instructive to see how the problem with trying to get rid of V2 can be seen in 
the fermion representation; in Appendix B the analysis will be done in considerable detail 
using the boson representation introduced in the next sub-section. Using the even-odd 
representation of Eg (|30|) , one can indeed rediagonalize the even-odd channels and absorb 
the Vi, V2 terms into Ho- This leaves the long time Green's functions of the even-odd 
channels unaffected but changes the short time behavior (see, for example, Nozieres and De 
Dominicis0). Thus 

(4(r)ce(0)) ^ {ct{r)cM) (36) 

for r small, of order the cutoff t^. Performing perturbation theory in V3 to second order we 
get a correction in U 

SU OC Vi {etc, ((c.C+) - (c+C„)) + C+Co ((ceC+) - (c+Ce)) } . (37) 

In general, (c^c) 7^ (cc"*") due to the nonlinear dispersion of fermions away from kp, in 
particular particle-hole asjTumetry. Thus, from Eq(^) we see that the asymmetry between 
the even-odd channels reappears in perturbation theory, due to the short time (high-energy) 
details. As a result we must retain the V2 term in the Hamiltonian. 

B. Bosonization 

In order to proceed it is necessary to find a representation that focuses on the essential 
low energy parts of the problem. Then, even if the problem is not exactly solvable, one 
can at least hope to be able to understand the physics and predict the low energy behav- 
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ior. The most commonly used representations are boson representations of the pseudo-one 
dimensional fermions. 

The basic strategy of bosonization is to try to mimic the low energy physics of the Fermi 
system, which can only be done exactly for a particularly simple system of exactly linear one- 
dimensional bands with a specific form of the cutoff. In more general situations, it is hoped 
(or, better, demonstrated!) that the high energy terms that are ignored — for example 
particle-hole excitations far from the Fermi surface — only serve to give finite renormal- 
izations of the basic parameters of the dominant low energy operators in the Hamiltonian. 
High energy properties — for example the fermion anticommutation relations — are, ipso 
facto, only handled approximately. What is more important, but unfortunately sometimes 
forgotten, is that terms that are formally irrelevant at low energies can, either on their own, 
when combined with other terms, or under canonical transformations, produce relevant or 
marginal terms that affect the physics. As we shall see, this is the case for the present 
problem. 

For now we will proceed in the conventional manner. Let us then start with the non- 
interacting Hamiltonian, Hg. If we are primarily interested in energies close to ep, we can 
extend the linear dispersion relation to all energies. Thus if we set the origin oi k at k — kp 
for convenience, we will get by Fourier transforming: 

/■oo 
dx 

where x is the conjugate variable to k and ^'j(a;) is the Fourier transform of Cjk- 

/oo f]U 
_^^e^^%-, J = 1,2 (39) 

Note that there are only right moving fermions in the system since there is only one Fermi 
point, i.e. one /c-value at which = ep. 

The operators at the impurity sites are given in terms of 

cj = J'^cj,^^j{0) J = 1,2 (40) 
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(38) 



which identifies cj as the x = creation operator in the one dimensional picture. Roughly 
speaking, "^jix) with x < corresponds to incoming s-wave electrons around site j while 
with X > 0, it corresponds to outgoing s-wave electrons. Thus the time reversal operator 
acting on "^jix) will give: 



(41) 



since it transforms incoming to outgoing electrons and vice-versa. 
At this point we can introduce the bosonic fields ^j{x) byEl: 



(42) 



with T^^ (X kp the characteristic cutoff frequency of the order the Fermi energy, and 



^j{x) = y/n (f)j{x) - / Ilj{x')dx' 

J —oo 



(43) 



with (pj and 11, satisfying appropriate commutation relations: 



{x),Ui{y)] = iSi. 



(44) 



\{T'^ + {x-yf) 

In the continuum limit Tc 0, this commutation relation approaches i5ij5{x — y). From 
Eq(^) and Eq(|42|) it can be seen that $j(x) transforms under time reversal as follows: 



f<|.,(x) = -<l>,(-x). 



(45) 



Expanding (f)j{x) and Ilj{x) in terms of their Fourier components. 



(j)j{x) 



Uj{x) 



°° de 



=° de\e\ 
°° 27T^/2\e 



(e)e'^" + 0;(e)e"'^" 



e 2 



_,0^.(e)e-- + z0+(e)e-^^" 



e 2 



and inserting these expressions in Eq(|43|), can be written as: 



°° de 



V2 



7re 



[e)e*^" + </)+(e)e- 



e 2 



(46) 



(47) 
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which involves only the positive energy parts. Subsequently inserting Eq(^7D in Eq 
then in Eq(B8|) the non-constant part of becomes 



j=l,2 ^' ^" 



e)e- 



and 



(48) 



with 



0,(e),0+(e')] = 5.,27r5(e - e') 
Finally using the standard expressionS 

1 d^j{x) 



(49) 



(50) 



the potential U is found to be: 



U 



27r 



9<l>i(0) ^ d^2{'^) 



dx 



dx 



■KTr, 



cos[$i(0) + $2(0)] + -^a 



27r 



a$2(o) d^i{o) 



dx 



dx 



(51) 



The form of U may be simplified by introducing Bose fields corresponding to collective modes 
for excitations that are even, $e(a^), and odd, $o(a;), about the center of symmetry of the 
pair of sites: 



$e= -^($l + $2) 
*o = ^ (*2 - $l) 



(52) 



The 0e,o and Ile,o can be defined equivalently. In terms of the new variables Ho remains 
in diagonal quadratic form (i.e. the indices j = 1,2 in Eq(^8D are replaced by j' = e,o), 
corresponding to free bosons and U becomes 



^ = m + cos V2$o + — =-(T^— 

V27r ox TTTc V27r ox 



(53) 

The symmetries of the problem are manifest in Eq(|53|), <l>e ^ "^'e + c corresponding to 
the global gauge invariance, Oz — o"z with $o — $o corresponding to the interchange 
symmetry, and $o — ^ '^'o + vr-\/2 corresponding to c- —Cj for j = 1, 2. Note that 
the even mode is completely decoupled from the impurity in Eq (|53D and will therefore not 



play a role. Formally it can be eliminated by a unitary transformation H AeHA~^ with 
Ag = exp[— ^^^^$6(0)]; we perform this transformation and henceforth only consider the 
potential Eq(|53D without the Vi term. 

In addition, Eq(p3D also has the symmetry: <l>e,o —>■ —^e,o and x —x corresponding 
to time reversal invariance from Eq(E3p. If the system were not time reversal invariant, 



then one could have V4 7^ 0. Indeed, the lowest order time reversal symmetry breaking 



term is sin 



such terms as gj^ 

ox 



v^$o(0) which is exactly the V4 term. However, we argued (above Eq(|^)) 
that V4 can always be chosen to be zero. But in the absence of time reversal invariance, 

V^$o(0) can also appear, essentially from nonlinear dispersion of 
the fermions away from kp and energy dependence of the scattering, that breaks the time 
reversal symmetry of the even or odd channels. Although these appear to be irrelevant, 
they cannot simply be eliminated because they generate a sin \/2^o{0) term after the 
unitary transformation of Eq(^) is performed. In order to eliminate such terms one has to 
pick a gauge or, equivalently a basis for the fermions (i.e. pick appropriate a, /3, 9 in Eq(|26|)) 
which creates a V4-term that exactly cancels the generated sin ^/2$o(0) term. In effect, 
one would thus obtain a set of almost time-reversal invariant low energy excitations, and 
our main results would still obtain. 

But danger lurks: even with full time reversal invariance similar terms to those discussed 
above will invalidate a related form of bosonization that we now discuss. It is tempting to 
find a way to get rid of the V2 term by a different choice of bosonization. One way to do 
this is to start with U in the form of Eq(pO|) and bosonize the fields Cek and Cok- In this case 
one has to introduce the fields \I/e(a;) and \I/o(a;) in an analogous way to \l/i(x) and ^1/2(2;) in 
Eq(ll): 

y/ZTCTc 

%{x) = ^e-^*^(^) (54) 

with $gQ(x) different from $e,o(a;) in Eq(^). Subsequently by introducing linear combina- 
tions of these fields 
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$6 = ^ - K) , (55) 



U would take the form: 



^ = m 7^ + ~7n + cos V2$b 56) 

V27r ox V27r dx ttTc 



with gauge symmetry under ^ + c; interchange symmetry under — o"^, $a ^ 

$a + — and time reversal sjTumetry under x —>■ —x, $a,b —^a,b- By then 

performing a unitary transformation using 



A' = exp 



,J_4„(0)+.-|-t,(0) 



(57) 



the first two terms in Eq (|56|) vanish so that the transformed Hamiltonian becomes 



A' HA'-' = Ho + —a, cos[V2$,(0)] + A^a,. (58) 



Noting that the second term is nothing but the third term in Eq(|53D, one might be tempted 
to conclude that the dynamics of the impurity is independent of V2. However, a term of 
the form cr^^^^ will appear in Eq(|5BD from a perturbation expansion in Aq and V2 which 
cannot be eliminated by the canonical transformation and is not irrelevant. Indeed, this term 
breaks the apparent symmetry of Eq (|58|) under $5 ^ — $5 which corresponds to $g 
and is not an exact symmetry in the presence of V2 or similar terms. In the Appendix B 
we consider a more general formulation of the two site problem which shows how, even if 
the part of the Hamiltonian that is symmetric under az — > —Cz is diagonalized fully before 
bosonizing, the energy dependence of scattering processes will, nevertheless, generally lead 
to terms which play a similar role to V2. 

We now analyze the effects of the even-odd symmetry breaking terms, such as V2, and 
other terms that are generated from these, by use of a Coulomb gas representation. As 
we shall see, in certain regimes, extra relevant operators are generated from the marginal 
operators that will delocalize the heavy particle. 
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III. COMPUTATIONS 



A. Path Integral Representation 

In this section, we analyze the bosonized Hamiltonian Eq(^) by a Coulomb gas repre- 
sentation and show how extra operators are generated which physically correspond to the 
impurity hopping together with one or more electrons. These will delocalize the impurity in 
regimes in which it was previously believed to be localized. (They will also generate extra 
important operators in Hamiltonians with " electron assisted tunnelling" like that analyzed 
by Vladar and Zawadowskii.) The Hamiltonian Eq(^3D is in a convenient form since it in- 
cludes, in a simple way, a term that breaks the artificial symmetry in the absence of V2, as 
well as a simple form of the impurity coupling to the Fermi sea. To treat Eq(|53|) we first 
perform a canonical transformation from H to H' using the unitary operator 



A = exp 

Then H becomes 



■j-..(0) + -|-....(0) 



(59) 



H' = AHA-^ = Ho + — cos[V2<l>o(0)] + {a+exp kV2Qo^oiO)\ + h.c] (60) 



where cr± = ^ {a^ ± iCy) and Qo = which as will be seen later is the effective charge 
for the hopping process. In Appendix A it will become clear why really is the charge 
transferred when the impurity hops between far away sites, by expressing Qo in terms of the 
scattering phase shifts. 

The Hamiltonian Eq(|60D is expressed entirely in terms of exponentials of boson operators 
which are particularly convenient for deriving a Coulomb gas representation. Note also that 
the even parts of the Bose field are completely decoupled from the odd parts and the impurity. 
Although the formally irrelevant operators in Eq(|53D have been ignored, their effects under 
the canonical transformation would only be to modify the coefficient of the V2 term, and 
to give operators that are still irrelevant and break no symmetries, although they would 
include coupling to the even part of the Bose field. We can however safely ignore these. 
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The correlations of the impurity position cr^, will not be affected by the canonical trans- 
formation as cr^ commutes with A. We therefore work with Eq(^) and drop the prime on 
H. Including the effects of spin a, we have, dropping the "o" (odd) subscript on $ 



Vn 



H = Ho + — J2 cos [72$'^ (0)] + A J a+ n exp \tV2Qo<^'' (0) 



TXT, 



c a 



h.c. 



(61) 



with 



~r Jo 2tc 



e e 



(62) 



and 



$"(x) 



de 



0-(e)e'- + 0-+(e)e- 



e 2 . 



-f3H 



(63) 
in the limit 



/o V27re 

We are interested in the zero temperature partition function Z = Tr 
(5 ^ oo. If we expand Z in V2 and Aq we obtain a sum over one dimensional "paths" from 
r = to /3. Each of these paths corresponds to a process in which the impurity hops between 
the sites at particular times shifting the phase of electron excitations, while at other times 
the electrons hop via the V2 term. Such a path is illustrated in Figure 2. For simplicity, we 
first work with a single spin species. Then we can write Z as 

00 00 2m 



n=Om=0 k=l \Cfe=±l/ 
dSi... dS2m / dri... dT2nZnm{{Ck}, {Si}, {Tj}) 

JO Jo Jo 



(64) 



with 

Znm{{Ck},{Si},{Tj}) = (0 



T 



. 



(65) 



where the product over different Cfe = corresponds to the two different terms in 



cos (-\/2$^ = I (^e*^* + e and the "fugacity" y of the electron hops is defined as 

V2 



y = 



2nTr 



(66) 
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We have taken the impurity to be on the "one" site at r = 0, this merely reduces Z by 
a multiphcative factor of two. Note that the signs of the ±i^/2Qo^ (r„) must alternate 
corresponding to the particle hopping back and forth, i.e. alternating cr+ and cr_ terms from 
Eq(^), and the sum is constrained to an even number of hops because the impurity begins 
and ends at the same site. In addition, we shall see that only terms with an even number 
of y "charges" will contribute, thus the sum J^l^i Cfc has to vanish. We have suppressed the 
dependence of $(x = 0) on the variable x and by $(r) we denote: 

$(r) = e+^°^$e-^°^ (67) 

with the expectation in Eq(^) taken with the ground state of Hq. The evaluation of the 
expectation value of the time ordered product in Eq(|65|) is particularly simple due to the 
bosonic character of $. First we observe that Eq(|67|) becomes 

de 



e""?. (68) 



/o VSvre 

Thus it is of the form $('r) = B{t) + B^{—t) where B{t) is a boson having commutation 
relations: 

[B (r) , B (r')] = 

and \b (r) , (-r') = / —e-'^^e-<^~^"^ = I{t - r') (69) 

for T > r', I{t — r') being formally divergent at small energies; but only the finite part 

/(r) ^ /(r) - /(O) = In (^) (70) 

will enter physical quantities. By the standard procedure of reordering the operators to bring 
B to the right and B^ to the left using the commutators in Eq(|69D, and noting that for 
zero temperature, only the ground state of Ho will appear on the right and left, we see that 
terms with the total charge J2kCk 7^ will give negative infinite terms in the exponentials 
and thus zero contribution to Z. Furthermore the partition function Znm can be written in 
terms of effective interaction between the charges with strength 
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W,,=2m,i{\r,-T,\) (71) 

where rj = Xj with qi = +Qo for the cr+ impurity hops 1 ^ 2; = Xj with qi = —Qo 
for the cr_ impurity hops 2^1; and Xj = Sj with q^ = Q for the electron hops, y. The 
partition function Znm is thus simply the Boltzmann factor for the charges interacting with 
the logarithmic potential Eq(|7TD, i.e. 

Z^m{{Ck},{si},{r,}) = e-'- (72) 

with 

2n 

m 

{k,k') \ '^c / 

2'^2n /U_T-|-l--r\ 

+2go E E (-1)' a in ' ^ ''^ (73) 

fe=l 1=1 \ ''"c / 

It is instructive to note that the exponential of minus the second term in Eq ([73|) has the 
same form as the amplitude (A^)^ in Eq(|^) which is nothing else than the square of an m 
particle Greens function. The reason it is the square of Am and not just Am is that here we 
have the product of the Greens functions of two sets of particles, the Ci^s and the c^'s. 

The partition function Z in Eq(0) with Znm from Eq([73D, is thus a 1-D Coulomb gas 
with logarithmic interactions between integer charges ±Cfc with fugacity y and "hopping" 
charges with fugacity Aq and charge ±Qo which must strictly alternate in sign. 

To take into consideration the effects of more than one spin species, we must modify 
Eq(^) to include the effective interactions between the impurity hops a± from each spin 
and the ±1 charges from the V2 term in H for each spin species. Since the Bose fields are 
independent, the interactions will be simply additive. Thus we must replace the "charges" qi 
in Eq(^) by vector charges ql with two components for spin-1/2 electrons. For the impurity 
hops we have 

q = Qo = (QcQo) 
or q = -Qo = {-Qo, -Qo) 

24 



for 1 — > 2 or 2 — s> 1 respectively, while for the electron hops, we have 

q= (CO) 
or q = (0,C) 

for spin t or J, electrons respectively with C, = ±1, depending on the direction of the electron 
hop. The total effective action now has the form 

Snm = Snm - 2m\ny - 2n In 

= -2^gl ■ g;; in f ^^~^'^^A - 2mlny - 2nlnAo (74) 

j>i \ J 

with the ordered times = Sj or Xj and the 2n charges Qi = ±Qo with fugacity Aq, strictly 
alternating in sign. The 2m charges in Eq(|7^ with fugacity y, need not alternate in sign 
but the sum of all their charges must be zero. 

B. Renormalization Group Analysis 

We are now in the position to analyze the behavior of the impurity in terms of the 
properties of the generalized Coulomb gas with the action of Eq(|7^. We are particularly 
interested here in whether or not the impurity can be localized. Thus we consider the effects 
of a small hopping rate Aq and analyze the Coulomb gas perturbatively in both Aq and y. 
Standard balancing of the energy and "entropy" terms indicate that y is exactly marginal 
— as it must be — , while the hopping Aq has renormalization group (RG) eigenvalue 

Ao = 1 - 2Ql (75) 

under rescaling of the small time cutoff, Tc, with the factor of 2 coming from the two 
spin species. Since for Qo > 1/V2, \o < 0, it would appear that the impurity hopping 
is irrelevant in this case, thereby leading to the conventional conclusion that a charge two 
impurity, which has Qo = (1,1) corresponding to the two spin channels, can be localized. As 
we shall see, however, extra charges will be generated under renormalization which invalidate 
this conclusion. 
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— * — * — * 

Since there are two types of charges ( and ±Qo, of which only the ±Qo are restricted to 
alternate in sign, there are various processes which can be regarded as composite charges. For 
example, a charge q — {Qo — 1, Qo) can be formed if an impurity hop and an electronic hole 
hop are close to each other. This process as well as the more general processes which generate 
charges {Qo ± n, Qo ± n) with n — 1,2,3,... were discussed physically in the Introduction. 
It is clear that while they do not exist in the original Hamiltonian they are generated under 
renormalization (or from perturbation theory). For example the hopping matrix element 
A(_i 0) of a process with charge {Qo — 1, Qo) will be generated under renormalization with 
magnitude proportional to yAg. Thus, in general, we must consider all possible composite 
charges and their effects on each other. 

We denote charges associated with general types of impurity hops: Q and —Q for a+ 
and (T_ hops, respectively, which occur at times r^; ( for purely electronic hops at times s^; 
and charge q for generic hops of either type at times with fugacity z. We need to analyze 
the effects of integrating out all pairs of charges with spacings between Tc, the cutoff and 
hence the minimum allowed spacing, and Tc{l + dl) with e' the time rescaling factor. Pairs of 
charges, qi and ^, which do not sum to charge zero, will generate composite charges + q2 
with fugacity ziZ2ki2TcSl with the kTcSl factor from the possible separations and ordering of 
qi and q2, with k = 2 unless both q's are impurity hops, in which case A; = if they are both 
(7+ and k — 1 (since only one ordering is possible), if one is (7+ and the other (T_. Neutral 
pairs, i.e. those with total charge zero, will not produce composite charges but will screen 
the interactions between the remaining charges. 

We can proceed as usual by considering the effects of one neutral pair on other charges, 

— * 

specifically on a charge q at time r. If the pair is ±C, i.e. purely electronic, then the effect 
of the two possible orderings cancel (up to modifying sub-logarithmic corrections to the 
interaction between remaining charges) and we thus ignore these. The interesting case is 
thus a pair of hops. The allowed orderings of a ±Q pair at times r ± Tc/2 with fugacities 
= A_^, depends on ^"^(t). Thus the interaction of this pair with charge q at time r is 
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6I,ir) = -2nQ ■ (76) 
T — r 

for |r — r| ^ r^, the appropriate limit for analyzing the renormalization of the long time 
interactions. Expanding in 51^{r) and integrating over the possible position, r, of the 
pair and the intra-pair spacing in the range Tc to Tc(1+(5Z) we see that there are contributions 
every time (7z{t) changes sign, i.e. at times Tj. This generates an effective interaction between 
q and all impurity hops, but not between q and purely electronic hops. For an impurity hop 

— * 

iQj at time Tj, the generated effective interaction is 

51.^^^ = ±4r,2A|4 • ? In \n - A 51. (77) 

— * 

This thus has the effect of modifying the interaction of each Qi with all other charges by a 
way that is equivalent to changing Qi by 

8Qi^-2Tl^lQ5l. (78) 

— * 

with zbQ the charges of the electron impurity impurity hop pair that have been integrated 
out. Since each Qi is of the form Qg plus an integer vector, we see that the net effect is just 
to change Qo by 

5Qo = -2t151 (Qo + n) A% (79) 

N 

with the sum running over all possible types Q of 0"+ charges, i.e. Q — Qo + N with N an 

integer vector; here and henceforth we use the abbreviated notation 

To this order in Aj^ and yj^, the fugacities for multi-electron hops which can have all integer 
vectors N except (0, 0), the RG flow equations are, after absorbing Tc's into A and y to make 
them dimensionless 
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dl 

dA 



N 



dl 

dl 



1 



N 

2 



Qo + N ) + 2 ^ y^, A^_^, 



2 



+ E Vn' Vn-N' + E ^N' ^N-N> (81) 
N'^Q N' 



with 



by the 1^2 interchange symmetry. As can be seen from Eq (|8lD all the multi-electron hop 
terms are irrelevant, thus we need only retain ?/(io) = y(oi) = y, from Eq(|8l|) we see that this 
is sufficient to generate all the composite charges with fugacities of order Aq times powers 
of y. 

From Eq(^) we see that, generically for two spin channels, there are at least three 
relevant operators for any Qo'. A([qj_ijq^]) and A([q^] ^gj^i) with [Qo] the fractional part of 
Qo are always relevant while either A([q^] [q^]) or A([q^]„i [q^]„i) or both will also be relevant. 
More generally, we arrive at the same conclusion as from the simple physical argument of the 
Introduction: in order to localize the impurity, more than four channels (including spin) are 
needed so that, if each channel, 7, is optimally coupled by a ^-integer Q^, then J2 > ^ 
and the impurity can be localized. 

The RG equations (Eq(^)) are quite different from those in the literature: if there were 
no V2 term, then the extra composite impurity hopping terms would not be generated, and 
the impurity could appear to be localized by just s-wave scattering. Note that the apparent 
Ci —>■ e*'^Ci symmetry when V2 = actually allows for charges Q = Qo + N with even integer 
vector A^, but these will not prevent localization. This point suggests that Qo should be 
defined up to an even integer (i.e. mod 2). In fact it will be seen in Appendix A that Qo is 
actually vr"^ times a phase shift which naturally leads to its consideration of mod 2. 

In the Hamiltonian considered by Vladar and Zawadowskii single electron assisted hop- 
ping terms that correspond to A(_i q) and A(o,-i) are present but their Hamiltonian has the 
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implicit symmetry c\ — ci, d\ —d\. Then only a subset of the can be generated, 
specifically those with + A^^^ odd ; again this artificial extra symmetry will change the 
behavior by limiting the number of possibly relevant operators. 

IV. DISCUSSION AND CONCLUSIONS 

In the previous section, we have seen that an impurity hopping between two symmet- 
rically placed sites cannot be localized unless it is coupled strongly to more than four spin 
and angular momentum "channels" . If the sites are nearby — as they must be if the bare 
hopping rate is to be appreciable — then one cannot use angular momentum channels, and 
must, instead, generalize the treatment. One way to do this, which shows directly the role of 
the irrelevant operators and relies on no symmetries other that the site interchange symme- 
try, is to use the one-electron eigenstates of the symmetrized electron Hamiltonian Hs which 
is the average of the Hamiltonians with the impurity on the two sites. The antisymmetric 
part then scatters the electrons between even and odd parity states of Hs- The analysis in 
this representation is carried out in Appendix B, with the same conclusions being reached 
as in Section III. 

The problems with most earlier treatments of the two-site system have been of two types: 
In many of the treatments, an extra ?7 (1) symmetry associated with the independence of 
the electrons which interact with the impurity at the two sites is implicitly assume d0. The 
V2 term that breaks this symmetry is marginal but it creates extra operators, particularly 
those which move one localized hole with the impurity, and these processes delocalize a 
charge two particle which had previously claimed to be localized if it interacts with only 
s-wave electrons at each impurity. 

Recently, there has also been a substantial literature on the relationship between an 
impurity with "electron assisted tunnelling" — i.e. hopping of the impurity concomitantly 
with the motion of one electron of either spin — and the two channel Kondo problem 
with the z-component of the Kondo "spin" being the impurity position and its x-component 



29 



the hopping, (i.e. our az and (Tx)- The two "channels" of the Kondo problem are then 
the two electron spin species which are exactly degenerate in the absence of an external 
field. The frequently used Hamiltonian0 for this problem was introduced by Vladar and 
Zawadowskii. However, they completely neglected a V2-like term. They derive the electron- 
assisted hopping term assuming that Qo ^ 1 and kpR ^ 1- Even in this case, using their 
numbers we find that the estimate of the amplitude of the electron-assisted hopping term 
they get is smaller by at least one to two orders of magnitude from the amplitude of A(_i q) 
that is generated from Eq(^) after renormalizing to / = 0(1): 

A(_i,o) ~ y^o ~ ^oT^PfV (83) 

with kpR < 1. In the case of kpR > 1, essential to get Qo = 0(1) which is the relevant 
situation for localization, their derivation of the electron-assisted hopping term breaks down. 
In our treatment and because of the existence of V2 we show how both A(_i^o) and other 
relevant terms (e.g. A(_i _i)) naturally arise. Thus, in any case, we believe we have here a 
more complete physical picture of the problem. 

As we have seen, these extra terms change the physics for large Qo and small hop- 
ping — the "weak coupling" (small J) limit in the Kondo language. Most of the recent 
literature@Jii"Ei has focused on the "strong coupling" behavior of the Kondo system i.e. the 
regime at low energies with parameters, particularly Qo, such that exchange is relevant and 
flows to large values corresponding to the particle hopping back and forth between the two 
sites. Novel non-Fermi liquid behavior has been found theoretically in this regime for which 
our considerations are not directly relevant. Nevertheless, the extra symmetry implicit in 
these treatments is potentially dangerous, indeed, as we will show elsewhere, terms that 
break this symmetry change the physics in the strong coupling delocalized regime, as they 
did in the weak hopping regime analyzed in this paper. 

One of the major motivations for the present work was the hope of gaining further 
understanding about the properties of a heavy particle that can hop on sites of a periodic 
lattice impeded (or in some regimes assisted!) by the coupling to a Fermi sea. This has 
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potential relevance for the mobility of muons in metalsB, the sharpness (or rounding) of 
X-ray edge singularities when the deep hole can move (albeit with a large bare mass), and 
possibly the properties of a heavy d- or f-electron band coupled to a light conduction band0. 
The main papers (e.g. reference 1 and references therein) on the behavior of a single particle 
in such a periodic system, suffer from some of the same problems as those on the two-site 
system: they treat a subset of the allowed operators and do not allow for the effects of others 
that may be generated. Not surprisingly, the conclusions of these papers are the same as 
for the two site case: that a charge two particle with only s-wave scattering can be localized 
while a charge one particle cannot be. In light of the present results, this conclusion should 
clearly be reexamined. 

In the spirit of the work of Sols and Guineai, one could treat each step of the particle 
motion — whether via nearest or further neighbor hopping — essentially independently and 
look at the renormalization of each such hopping term separately, including the possible 
motion of electrons with the particle, by the methods outlined in this paper. From this 
approach the following conclusion would be immediate: that the particle cannot be localized 
unless all the hopping processes are irrelevant, and this can only occur if there are more 
than four electron channels involved for every possible hopping process and the coupling 
is sufficiently strong. Thus with just s-wave scattering, a particle cannot be localized, in 
contrast to the conclusion in the literature0i. Instead, the particle will move around with a 
screening cloud of electrons in tow. 

Unfortunately, there are problems in extending the two site results in this way. The 
primary one is the spatial structure of the system and the lack of independence between the 
electrons involved in hopping between different pairs of sites. For any finite number of sites, 
the electrons can be treated as essentially one dimensional at sufficiently low energies, and 
generalization of the present methods can be used to categorize all the operators. We have 
explicitly carried out this procedure for a simple case of three sites symmetrically arranged 
in a triangle and find that the same conclusions are obtained as in the two site case. We 
believe that this should likewise hold for any finite number of sites. However for an infinite 
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lattice of sites, the electrons must be considered to be fully three dimensional and the limit 
of small bare hopping rate that we have studied perturbatively, may not be exchangeable 
with the limit of an infinite number of sites. In particular, one has to worry about at least 
two effects. First, even if the particle is moving very slowly, there will always be particle- hole 
excitations with group velocities slower than the particle; these can perhaps not be treated 
in the same manner as the rapidly moving excitations. Second, one could argue that there 
are an infinite number of scattering "channels" involved because of electrons near each of the 
sites. Perhaps these might make it easier to localize the particle, although it is not clear how 
this could come about. Conversely, they might somehow interfere and prevent localization 
even if there are many scattering channels at each site. We must, unfortunately, leave these 
questions for future study. But we should note, that if the preliminary conclusion about 
the difficulty of localizing a heavy particle is correct, it may have implications for some of 
the suggested possibilities for interesting new physics with a heavy f-electron band coupled 
to a light conduction electron band@. Again, we leave this, as well as the possibility of 
interesting effects on the electrons caused by a delocalized impurity, for future work. 

The present analysis, although not introducing new techniques, has, we hope, made the 
physical picture behind the competition between orthogonality catastrophe and hopping 
clearer and simpler. In addition, the pitfalls of the standard combination of bosonization 
techniques and "large" canonical transformations have been brought out and should be 
heeded by workers on other problems in this area. 

ACKNOWLEDGMENTS 

We thank Bert Halperin, Andy Millis and Igor Smolyarenko for useful discussions. This 
work was supported in part by the NSF via grant DMR 9106237. 



32 



APPENDIX A: 



In this Appendix it is illustrated how one can get a Coulomb gas representation of the 
partition function Z directly from the original form of the Hamiltonian Eq(|3^), without 
bosonizing. At the same time the physical interpretation of Qo as the transferred electron 
screening charge instead of simply being the interacting potential V3 — which is the naive 
result of bosonization — will become apparent. The method used here has been applied 
to systems very similar to ours in the pasti'Bil so we will not go into the details of the 
calculations. 

It is straightforward to see that using the Hamiltonian in Eq (^2|) one can get an expression 
for the partition function Z similar to Eq(|64|): 



n=0 m=0 

where ?/ oc V2 and 



'Z'nm ({-Si}, {'^0'}) 



ds\ . . . 

JQ 



dS2m / dTi... 

JO Jo 



dT2nZnm {{Si} , {t,}) (Al) 



T 



cl{s{)C2{Si) . . . 4is2m)Cl{s2m) CXp f - ^ dr' H'^ (t') 



In EqM) by H'^ we denote 



. (A2) 



H'^{t) = H,+ {V, + Vs) cfc, + {V, - Vs) C+C2 + 2\/3^(r) (c+C2 - c+c. 



(A3) 



with 



9{t) 



for Te{T2k,T2k 



+1) 



A; = 0, 1, . . . , n 



1 for Te{T2k-i,T2k) k = l,...,n 



(A4) 



Since the V2 term contains both C2C1 and C1C2 which can be generated at any time Si, there 
is no constraint in the order in which they appear in Eq([A^) . However, due to the number- 
conserving character of H'^{t) there is a constraint of having, for each z, as many cf as Cj's 
in the time ordered product of Znm- For convenience we pick |0 > to be the ground state of 
Ho + iVi + V3) CiCi + {Vi — V3) c^C2. We can see that Znm represents a Green's function of m 
particles, types "1" and "2" , in a time dependent potential which changes its value when the 
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impurity hops at times Tj. This, in the language of the Kondo proble m§, is equivalent to a 
path in which the impurity spin flips at times while the electrons flip their "spin" (i.e. 1-2 
index) at times sj. The amplitude for such a process is a product of two independent parts, 
i.e. the amplitude of the 1-particles and that of the 2-particles. Now since the "interaction" 
term V3 is really only a one particle operator, all the diagrams will be either closed loops 
or else will end at times Sj, so we can thus treat them independent 1}0. The closed loop 
contribution is simply Zno which is identical to that of the Kondo problem. 

In order to get the other contribution let us first consider m=l. Then the amplitude corre- 
sponds to a one-particle Green's function in the presence of the time dependent Hamiltonian 
which can be shown to have the following long time (i.e. all time differences |r — r'| ^ r^) 
solutionill: 



1 " 

G{s,s') oc — - n 

k=i 



ij2k - S) {T2k-l - S' 



{r2k - s') (r2fc-i - s) 



(A5) 



with the phase shiftE 



6 = - arctan + nQ{-f3) (A6) 



with the Heavyside step function and 

/3 = 1 - V^i tan + ^ (1 + tan^ o) (v^ - V^) (A7) 

with 6 related to the short time behavior of Go (the propagator of Ho) including the existence 
of bound states and particle-hole asymmetry0. Note that the extra 0-function, which in 
the past had often been ignoredi'0 and was only recently introduced explicitlyil makes S 



defined between — | and |. 



From Eq (|A6|) it becomes clear that although Vi may be formally decoupled from the 
impurity in the bosonized version of the problem (see Eq(^)) it does renormalize the expo- 
nent 6. This leads us to the conclusion that indeed 5 is a non-universal quantity — not in 
general simply related to the phase shifts for scattering off a static impurityi. 

More generally, since electron operators anti-commute, the amplitude for m particles is 
a determinant of G{si, s'j) terms0 (see Eq(^): detG (^Si, s'j'j, where Sj and s'j are times at 
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which the i-th 1-particle was created and the j-th 1-particle was annihilated, respectively. 
Expanding the determinant (over i,j) out, we get 



detG'(s„s;.) 



nn 

i k=l 



[r2k - Si) [T2k-1 



[T2k - Si) {T2k-1 



(A8) 



We can now bring this in a form closer to Eq(^) if we introduce a number Q which is +1 
for Si and -1 for s[. Then, quoting the result for the closed loop amplitudeBS 



exp 



,1+1' 



In Ir^ - Ti' 



(A9) 



and using Eg (|A8|) , after taking the square of all the above amplitudes to take into account 
the two channels (1,2), Znm becomes: 

2 2n 



Znm {{Ck}, {Si}, {Tj}) = exp 



\ " / («,/') 



In |r/ - r// 



X exp 



X exp 



2m 



+2 CkCk' In |sfc - Sk' 

{k,k') 
^ 2m 2n 

IT 



(AlO) 



This is equivalent to Eq(^ 



k=l 1=1 

if ^ ^ Qo- The essential reason for the appearance of an 
arctangent of the potential in Eq([A6|) is because the singular effects of higher order terms 
in V3 are included. Eq([A^) makes 6 finite even if V3 is infinite. Finally, from Friedel's sum 
rule, it would seem that ^ corresponds to the electronic screening charge moved when the 
impurity hops from site to site, thus Qo would be just this screening charge. One should 
emphasize, however that this correct only in the limit of large site separation for which V2 
is small and 6 is not appreciably renormalized from the phase shift off a static impurity at 
one of the two sites. 

In conclusion, we have seen here that \Qo\ < |, since — f < 5 < ^. However this is 
essentially equivalent with the results of this paper in which we have shown that for any 
initial bare value of Qo the most relevant process that will dominate in delocalising the 
impurity, is the one that has charge q with absolute value 
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|g|=min{|[QJ|,|[Q„]-n|}<l (All) 

with [Qo] the fractional part of Qo (see discussion after Eq(pTD). The physical picture 
presented in this paper provides an intuitive way of understanding why the branch of 8 
implied by Eg (|A6|) is dominant. This was not apparent in some of the earlier wor 
and is the source of some of the erroneous claims about particle localization^'!. 



APPENDIX B: 

In this Appendix, we show how the site problem can be analyzed generally in terms of 
"channels" even in the absence of any symmetry except the equivalence of the two sites. 
In addition, we will see how bosonization in a representation in which no V2-like electron 
hopping term can exist will still, if handled carefully, yield terms which play the same role. 
Furthermore, the possibility of exchange of "charge" between various channels will be found, 
yielding another, albeit related way, that charges can be reduced, and localization impeded. 
We consider the electronic Hamiltonians with the impurity at either site one or two, denoted 
ll\ and if 2, and diagonalize exactly the symmetrized Hamiltonian 

iis = \ {Hi + H2) . (Bl) 

The full Hamiltonian is then written as 

H = Hs + a^HA + A.a, (B2) 

with the antisymmetric part 

Ha = \ {H, - H2) . (B3) 

With fermion operators that diagonalize Hs, we see that neither V2- nor Vi-like terms can 
occur. At each energy, there will be a countable degenerate set of scattering eigenfunctions 
of Hs, which we label by A; = e — ep, setting vp = 1] and for the states that are even (e) under 
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interchange 1 ^ 2 we use a "channel" index rj which we can choose later for convenience, 
while we use an index 9 for the odd (o) states. Thus 

Hs = J2J2 ^4neCkne + kcteo^keo (B4) 

k V k e 

and the antisymmetric part has the general form 

Ha = Y1 r,7,e(A;, k') ctneCk'eo + h.c. . (B5) 

k,k' 

Since we are interested in the behavior near the Fermi surface, we can choose linear combi- 
nations of the even states at A; = and likewise the k = odd states, so that 

r,,e(0,0) = 5^e7, (B6) 

is diagonal in rj, 6; then in this representation we can denote 6 also by t]. We now form the 
operators 

^r,e(o) = r— (^rike{o) (B7) 

Y ^e(o) k 

with appropriate normalization factors Ai'e(o) as in Eq(p3D, and define 

= ^ (^,e - *,o) . (B8) 

The Hamiltonian then takes the form 

H = Hs + a,Y.l^ (^ti^^i - ^-12^.2) + (J.M (B9) 
»? 

where the correction term M. involves the deviations of F from the diagonal form for k,k' ^ 0: 
M=Y.Y1 (rw'(A;, k') - 7^(5^^0 (4^eCfc'^'o + h.c) . (BIO) 

If we bosonize the ^1^,71,2 following the prescription used in Section II. B, then in the absence 
of M., the Hamiltonian just takes the simple form 
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with 

7n d^,o (0) 

i7, = ir, + -^a.^^ (B12) 

where is the kinetic energy of the rj bosons. We then have for each channel rj an 
independent charge 

Q,o = ^, (B13) 

71 

SO thatii 



a 



' TjO 



In the absence of M., the particle would thus be localized if fto > 1. 

Analysis of the form of M., shows, with the r] basis chosen to vary slowly with k, the 
existence of formally irrelevant terms like 

5<f,o(0) 



■ no \ , 

(Tz i COS 

OX 



v^$,,o(0) , (B15) 



the ^ essentially arising from terms in Jvl linear in k and k' . Under the canonical transfor- 



mation that eliminates the cr^ — terms in Eq (pi2| ), these will generate cos a/2$^o (0) 
terms which are of exactly the same form as these that would have arisen from a V2 term 
originally. These terms create integer charges which can then combine with the charges 
to give effective charges of Qr? — with integer n^. Because of the choice of the 77 fermions, 
these terms no longer have quite the interpretation of moving electrons with the impurity. 
Physically, this is quite simple: in the correct basis, the electrons do not need to be moved, 
they will do so on their own due to the change in the electronic part of H as the impurity 
moves. 

In addition to V2-like terms, terms of the form 

e^;^(*--*Ve) sin sin (^"^ (B16) 

will also be generated under the canonical transformation from F^^/ terms linear in k. These 
will create | charges in the formerly-decoupled even channels. For any pair 77, 77', it can be 

38 



seen that the effective charge squared of combining this process with a hop has a contribution 
to CKo from these two pairs of even and odd channels of 

/ 1\2 1 / 1 

[Qv-^) +4- (BIT) 

Since the resulting contribution from each even-channel charge is |, it can be seen that the 
effective charge squared of each channel cannot, in the general case, be reduced below |, i.e. 
the same result as in the absence of the channel mixing terms. Therefore these terms do not 
change the conclusion of the earlier discussion. Nevertheless, the presence of channel mixing 
terms will complicate the analysis of the many site problem and could perhaps change the 
physics in a spatially extended system. 
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